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21 Abstract 

22 

23 We compare estimates for the ion fluxes of twelve expected constituents of the lunar 

24 exosphere with estimates for the ion fluxes ejected from the lunar stuface by solar wind 

25 ions and electrons. Our estimates demonstrate that measurements of lunar ions will help 

26 constrain the abundances of many undetected species in the Itmar exosphere, particularly 

27 A1 and Si, because the expected ion flux levels from the exosphere exceed those from the 
23 surface. To correctly infer the relative abundances of exospheric ions and neutrals from 
2 ? Kaguya Ion Mass Analyzer (IMA) measurements, we must take into account the velocity 

30 distributions of local ions. The predicted spectrum underestimates the measured levels of 

31 relative to other Ixmar ion species, a result that may suggest contributions by 

32 molecular ions to the measured rates. 

33 
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1. Introduction 


35 

36 The Earth’s Moon is surrounded by a thin exosphere many of whose expected 

37 neutral constituents remain undetected [Stem, 1999]. The unknown constituents of the 

38 lunar exosphere can be probed by measurements of pickup ions, which are often easier to 

39 detect owing to the high sensitivity of ion mass spectrometers. Many ion species have 

40 been detected at the Moon although their neutrals have not. Lunar ions have been 

41 observed by the AMPTE [Hilchenbach et al, 1991], WIND [Mall et al, 1998], and 

42 SELENE (Kaguya) [Yokota et al., 2009; Tanaka et al., 2009] spacecraft. Such ions can 

43 be created not only by photoionization of exospheric neutrals, but also directly sputtered 

44 from the lunar surface during its bombardment by solar wind ions and electrons [Blpbic 

45 et al, 1991; McLain et al, 2011]. 

46 The relative importance of these two soiu-ces of ions, the surface and the neutral 

47 exosphere, is imclear. In principle, ions from the surface and ions from the exosphere can 

48 be distinguished by their distinctive energy spectra: the surface ions are nearly 

49 monoenergetic, while some exo-ion species may have a more extended energy range 

50 [Yokota and Saito, 2005], depending on the ratio of ion gyroradius to neutral scale height 

51 [Hartle and Killen, 2006] .However, it is difficult to quantify the different ion sources 

52 directly from measurements because surface ions can only be detected by spacecraft on 

53 the side of the Moon that lies in the direction of the solar wind electric field, whereas ions 

54 from the exosphere are more widely distributed [Hartle and Killen, 2006\. Therefore, 

55 self-consistent models for the expected neutral and exo-ion abundances are required to 
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56 clarify the likelihood of inferring neutral abundances for each species from ion 

57 measurements. 

58 Abundances of twelve ion species of exospheric origin (He"^, C^, O"^, Na^, Mg"^, Si"^, 

59 A^, S^, Ca"^, K*, Ti"^ and Fe^ are estimated herein, most of them for the first time. Our 

60 predictions are based on models for lunar exospheric neutrals [Hartle and Thomas, 1974; 

61 Sarantos et al., 2010; Sarantos et al. 2012] that were validated with existing 

62 observations (detections for He, Na and K; upper limits for other species). With these 

63 production rates, dependent on altitude and zenith angle, we initialized an ion transport 

64 model [Hartle et al, 2011] to quantify fluxes of lunar pickup ions from exospheric 

65 sources. These rates were compared to the expected rates of secondary ions emitted from 

66 the surface by sputtering and electron stimulated desoiption. Last, because ions near the 

67 Moon may be highly directional, we studied how exo-ions map into a detector that does 

68 not have a 4 ti field-of-view. This last part is specifically relevant to interpreting published 

69 measurements by Kaguya, because its time-of-flight Ion Mass Analyzer (IMA) is turned 

70 towards the Moon with a 27t field-of-view [Yokota et al, 2009; Tanaka et al, 2009; Saito 

71 et al, 2010]. 

72 

73 2. Observations 

74 Lunar ions have been observed far downstream of the Moon and, more recently, 

75 near the Moon. The first ions of limar origin were detected by the time-of-flight 

76 spectrometer onboard the AMPTE spacecraft; they were dominated by peaks 

77 corresponding to Si"^ and/or A^, with a secondary peak at [Hilchenbach et al, 1991]. 

78 Lunar ions were also detected during lunar swingbys of the WIND spacecraft well 
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79 upstream from Earth, yielding once again clear peaks due to O'*', sr and Ar, with oxygen 

80 being most prominent {Mali et al., 199S]. The observed could be attributed neither to 

81 magnetospheric nor to interstellar origin. More recently, pickup ions were observed at 

82 altitudes ~100 km above the Moon from the MAP-PACE-IMA spectrometer onboard the 

83 SELENE (Kaguya) limar orbiter. When the Moon was in the solar wind, Yokota et al. 

84 [2009] reported the detection of He'*', C^, O^, Na'*', and K'*'. In addition, when the Moon 

85 was in the magnetosphere, Tanaka et al. [2009] reported the presence of K*', He^, He'*', 

86 C'', O'*', Na^, and Ar'*'. Such ions could be a mix of solar wind, outflow from Earth’s 

87 ionosphere, and of lunar origin. 

88 Ions desorbed from the lunar surface can be constrained by laboratory experiments. 

89 When lunar soil simulants were bombarded by solar wind-like ions in Secondary Ion 

90 Mass Spectrometry (SIMS) experiments, significant fluxes of secondary lunar ions were 

91 ejected (-10^-10^ ions cm‘^ s'* if the typical solar wind flux of 4x10® ions cm'^ s'* is 

92 assumed) [Elphic et al., 1991]. Yields were higher for elements with lower ionization 

93 potentials such as K'*', Na'*', Ca'*', and aE. Experiments of electron impact onto Na and K 

94 bearing glasses indicated that Electron-Stimulated Desorption (ESD) produces ions 

95 {McLain et al., 2011]. Measured ESD yields for 1 keV electrons were ~ 10'® 

96 ions/electron, so the typical fluxes of electrons in the solar wind could result in a flux of 

97 ~10 total ions (summed over all species) per cm' s' from the lunar surface. Hence, the 

98 expected fluxes of surface ions from ESD and sputtering may be conqiarable in 

99 magnitude. ESD yields were found by McLain et al. [2011] to depend on substrate 

100 temperature, and hence on solar zenith angle. Owing to flieir dependence on the solar 
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101 wind incidence angle, sputtered ions will also have a solar zenith angle dependence, as 

102 will ions originating from the exosphere. 

103 

104 3. Integrated model of lunar neutrals and ions 

105 The lunar atmosphere can be produced by a variety of source processes; thermal 

106 and photon-stimulated desorption (PSD), micrometeroid impact vaporization, and solar 
lO”? wind sputtering {Stem, 1999]. The resulting gas populations exhibit a mix of thermal (or 

108 accommodated to the lunar surface temperature) and non-thermal velocity distributions. 

109 The altitude dependence of exospheric neutrals, a reflection of the ejection processes, is 

110 expected theoretically to control the distributions of resulting ions above the Moon 

1 1 1 [Hartle and Killen 2006; Hartle et al., 201 1]. 

112 The exospheric abundances of neutral He, Na, and K are known from 

113 measurements. We do not know the abundances of other expected constituents of the 

114 lunar exosphere but upper limits exist for many [Stem, 1999]. Wur 2 et al. [2007] and 

115 Sarantos et al. [2012] estimated abundances for these undetected species by assiuning 

116 thermal desorption, impact vaporization, and sputtering sources, respectively. Both of 

117 these studies concluded that models for many refractory gases are not well-constrained by 

118 existing measurements. In their study, Sarantos et al. [2012] found that; a) the published 

119 upper limits for Mg and Fe greatly exceed the expected production by micrometeoroids 

120 and the solar wind; b) the modeled Al and O column abundances from these sources are 

121 conq)axable to present upper limits; c) the modeled column abundance for Si and Ti is 

122 higher than published limits by a factor of four to six, and d) the calculated Ca abundance 
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123 


was much higher than the upper limit. The low upper limit for Si and Ti was due to 

124 measurements of lines corresponding to excited states that are unlikely to be populated. 

125 Additionally, the low upper limit for Ca could be interpreted as consistent with increased 

126 loss to condensation and molecular formation in the impact cloud. Because of significant 

127 uncertainties in microphysical parameters for the sources and losses, Sarantos et al. 

128 [2012] concluded that the predicted abundances of exospheric refractories are probably 

129 upper limits. 

130 Figure 1 presents the estimated distribution of twelve limar neutral species. 

131 Estimates for He were provided by models of Hartle and Thomas [1974], and estimates 

132 for Na by Sarantos et al. [2010]. Models for these two species have been validated with 

133 measurements. Estimates for the gaseous abundances of the main regolith constituents 

134 (O, Si, Al, Ca, Mg, K, Ti, and Fe) were obtained by Sarantos et al. [2012]. This model 

135 assumed that the siuface reservoir is uniform, that K is released by PSD with cross- 

136 sections and temperatures similar to that of Na, that refractories are released by the 

137 combination of ion sputtering and impact vaporization, and that all particles stick upon 

138 contact with the surface. Models for those species compared favorably to upper limits 

139 except for Ca; for that species we produced ions consistent with the upper limits [Stem, 

140 1999]. Because we do not have appropriate models for the production and loss of some 

141 expected volatiles, such as C and S, we used a spatially uniform Chamberlain [1963] 

142 model assuming the ejecta temperature of 400 K and the exospheric abundance at low 

143 altitudes given by Stem [1 999]. We do not present a model for Ar. 
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Pickup ions formed from the ionized neutral exosphere were tracked in the 

145 interplanetary electric and magnetic fields. Photoionization rates for these species were 

146 estimated from cross-sections by Vemer et ah [1996]. The motion of lunar ions under 

147 ExB drift can be obtained analytically (i.e., using constants of motion) when they travel 

148 in spatially uniform magnetic and electric fields. Under these conditions, Hartle et al. 

149 [201 1] employed the Vlasov equation with an ion source term to derive an expression for 

150 the phase space density, / of pickup ions that accounts for the three-dimensional 

151 distribution of the neutral source gas. With this approach, a given ion velocity vector, v, 

152 arrives at a location r with a probability density y(r,v) that relates to the neutral density 

153 back to its point of origin, to, through the expression: 

154 £ V(r„)(l) 

' d allCwhids 

(r.v)^(r,,v,=0) 

155 Here rg is the ion gyroradius; R the photoionization rate; = [E x B| / the 

156 bulk drift velocity of an ion in magnetic and electric fields B and E, respectively; N(re) 

157 the neutral density at each point of origin; Vy the instantaneous ion velocity along E at 

158 location r, and V;^ is the instantaneous ion velocity along the drift direction. The delta 

159 function indicates that /is a function of just one component of the ion velocity because 

160 under these assumptions ions move perpendicular to the magnetic field (vz = 0) and 

161 remain in a ring velocity distribution. Some velocities in Eq. (1) are less probable to be 

162 observed because they originated farther upstream, where the neutral density and the 

163 production rate of new ions are small. The summation includes all cycloids that 

164 originated several gyroradii upstream and arrive at the point of interest with the same 

165 velocity, v. In practice, only one to two cycloids must be included in the summation 
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166 because the gyroradii of lunar ions greatly exceed the neutral scale heights. Eq. (1) with 

167 N from our neutral models can be used to numerically compute how ions from the lunar 

1 68 exosphere map into an ion spectrometer. Integrals of the distribution function (1) over the 

169 allowed velocity space yield densities, velocities, and macroscopic fluxes at any location 

170 above the surface for arbitrary Interplanetary Magnetic Field (IMF) and solar wind 

171 directions. 

172 The formulation adopted here (Eq. 1) assiunes that the solar wind flow direction 

173 and magnetic field are spatially uniform, assumptions that are good approximations in the 

174 lunar dayside, but not in the lunar wake nor near the highly electrically charged lunar 

175 terminators. The simplifying assumptions of the ion transport model ignore the effects of 

176 wave-particle interactions, which may scatter ions in velocity space and remove them 

177 from a ring distribution, as well as the presence aroimd the Moon of magnetic anomalies. 

178 Given these limitations, our estimates are valid only for the limar dayside. 

179 

180 4. Model Predictions 

181 Figure 2a presents estimates for the ion flux of 12 expected exospheric constituents. 

182 These estimates were produced at 100 km, the orbital altitude of Kaguya, for a Parker 

183 spiral IMF configuration (field in the ecliptic, 45° off the Sun-Moon line), with a 

184 component away from the Sun and a magnitude of 8 nT. The assumed solar wind velocity 

185 was 400 km/s in these simulations. A realistic range for the expected fluxes of exo-ions 

186 can be obtained if estimates of omni-directional flux at two locations above the Moon are 

187 provided: at the subsolar point (blue), and at a point above the surface where the flux is 

188 maximum (red), a location away from the subsolar point in the direction of the electric 
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189 field (see Figure 3). Fluxes fi-om the exosphere are compared in Fig. 2b to the secondary 

190 ion flux leaving the surface [Elphic et al., 1991] due to solar wind bombardment at a flux 

191 typical of the solar wind at 1 AU, 4 x 10* cm'^ s'\ Based on these predictions, we 

192 conclude that the exo-ion rates for many of these species, especially and Si"^, greatly 

193 exceed the secondary ion rates. Hence, it should be possible to constrain the imknown 

194 levels of some lunar neutrals with ion spectrometer data. 

195 Now consider the spatial distribution of ions originating from the exosphere and 

196 their arrival into an instrument. Here we focus on the two best-specified neutral species in 

197 our model. He and Na. The exo-ion flux at an altitude of 100 km above the limar dayside 

198 is plotted in Figure 3. Shown on the left (panels 3a and 3c) is the omni-directional flux of 

199 exospheric He^ and Na^ as a function of Selenocentric Solar Ecliptic (SSE) longitude and 

200 latitude, with the subsolar point denoted by a white dot. Fluxes of dayside He^ peak near 

201 the terminator because He neutrals are more abundant on the cooler terminator [Hartle 

202 and Thomas, 1974]. The flxix of exo-ions peaks in the quadrant that lies in the direction 

203 of the electric field because at these locations we probe the lower volume (<100 km) of 

204 the exosphere (see Fig. 4 of Hartle and Killen [2006] for a visual explanation). 

205 The actual flux calculations are contrasted in Fig. 3b and 3d to the flux of the 

206 same two species that would be observed by a 2 ti detector pointing nadir - like the MAP- 

207 PACE-IMA spectrometer onboard Kaguya. Ions firom the exosphere map inefficiently 

208 into the nadir-pointing mass spectrometer in file quadrant against the direction of the 

209 solar wind electric field. This is because fi-eshly picked-up ions lie near the electric field 

210 direction, which for this IMF configuration is excluded fi'om the instrument’s field of 
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21 1 view in the region over soutiiem latitudes. A higher fraction of the total exo-ion flux is 

212 recorded elsewhere, with the portion of phase space density seen varying as a function of 

213 spacecraft location. The region northward of the white dotted lines in Panels a-d is the 

214 region over which both surface and exospheric ions map into the instrument. Counts 

215 below this line can be uniquely identified with an exospheric origin. Integrating the 

216 expected signal from exospheric and siuface ions over the Sun-Moon meridian to 

217 approximate the geometry of a polar orbiter like Kaguya, the model estimates that 

218 exosphere > 7 X Na‘^i„,/oce and Af exosphere > 30 X Al'*' surface for the ElphiC 6t ol [1991] 

219 yields and a solar zenith angle dependence for sputtered ions. Therefore, both local 

220 (instantaneous) spectra and spectra integrated over a dayside pass by the Kaguya IMA 

221 can be used to determine the relative composition of the exosphere. 

222 Our simulations indicate that some constituents of the exo-ionosphere 

223 map into Kaguya IMA better th an others, biasing the spectrum measured for a given orbit 

224 and IMF configuration. For example. Figure 4a shows the actual (omni-directional) 

225 Na'^/He^ flux ratio as a function of selenocentric location at an altitude of 100 km. The 

226 expected Na^ flux for a Parker IMF is about twice that of He"^ at subsolar locations, while 

227 the He"^ flux exceeds that for Na"^ by about a factor of ten on some locations near the 

228 terminator. Figure 4b shows the relative Na^/He"^ flux ratio for a virtual 2 ti nadir 

229 instrument like Kaguya IMA. Over large portions of the hemisphere on the side of the 

230 electric field, IMA would report a correct ratio, but, near the subsolar point, the relative 

231 abundance of Na'^/He'*^ measured by IMA would be higher than actual. Panels c and d of 

232 Figure 4 clarify that whereas the velocity distribution of Na"^ is sharply peaked, that of 

233 He"^ is wider, and more of it is excluded from the instrument over the subsolar point. In 
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234 contrast, over other locations, where Kaguya IMA cannot observe near the direction of 

235 the E-field, it would record a relatively higher portion of He"^. Because of these effects, if 

236 spectra integrated over a whole meridian are studied (Fig. 4e), the amount of He^ relative 

237 to Na"^ would be overestimated on the dayside for this IMF configuration. Therefore, 

238 Kaguya IMA spectra must be carefully interpreted orbit by orbit before accurate 

239 estimates of relative abundance can be made. 

240 5. Conclusions 

241 

242 Both siuTace and exospheric ions contribute to the lunar environment. Tr, Fe^, 

243 Mg"^, and especially Ca"^ are mainly ejected from the surface. For the other species studied 

244 here, ionization of the exospheric constituents could result in fluxes that significantly 

245 exceed the surface production rate. In our model, which predicts levels for neutral A1 and 

246 Si near the upper limits observed by telescopes, Al"^ and Si^ would be the dominant 

247 pickup ions species around the Moon, similar to what was observed by the AMPTE 

248 spacecraft [Hilchenbach et al, 1991]. Hence, we predict that A1 and Si should be the most 

249 easily constrainable neutral species by ion measurements. 

250 We cannot explain the reported overabimdance of observed relative to other 

251 species (roughly 3:1 count ratio between peaks corresponding to masses 16 and 28 in 

252 Mall et al, 1998) with neutral oxygen release at levels expected from impacts and solar 

253 wind sputtering. This conclusion is reached even though our neutral model prediction for 

254 O is essentially at the upper limits from measurements by Apollo 17 of the oxygen triplet 

255 [Table 2 of Sarantos et al, 2012 and references therein]. No was detected in sputtered 

256 ion experiments [Elphic et al., 1991] and little was detected in ESD experiments [McLain 
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257 et al, 2011], making a surface origin for unlikely. If a lunar origin of the observed 

258 is posited, volatile molecular ions such as OHT, H 2 O'*’, CO* or CO 2 * may be contributing 

259 to the O* and C* measured by Kaguya and WIND, a possibility suggested by Tanaka et 

260 al. [2009] that would help mitigate the model-data discrepancy for O*. 

261 Our estimate of lunar O* fluxes, ~ 4 - 6 ^ 10^ ions cm'^ s'*, suggests the 

262 possibility of significant contributions by Earth outflow to measurements obtained with 

263 the Moon inside the terrestrial magnetosphere. A model of die O* outflow from Earth 

264 [Glocer et al., 2009], run for conditions prevailing on April 20, 2008, the day of the 

265 observations reported by Tanaka et al [2009], predicts that ~3 x 10^ ions cm‘^ s"’ reach 

266 the Earth’s magnetic lobes at the distance of 60 Earth radii. Earth’s polar wind ions may 

267 affect the lunar environment even more significantly during geomagnetically active 

268 conditions. 

269 Discrepancies found between the simulated and observed ion spectra can further 

270 our understanding of the lunar environment, especially as new measurements of limar 

271 ions by die Acceleration, Recormection, Turbulence, and Electrodynamics of the Moon’s 

272 Interaction with the Sun (ARTEMIS) [Sibeck et al., 2011], and of hmar neutrals by the 

273 Lunar Reconnaissance Orbiter (LRO) and the Lunar Atmosphere and Dust Environment , 

274 Explorer (LADEE), become available over the next couple of years. 

275 
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337 Figure Captions 

338 Figure 1. Neutral profiles above the subsolar point that are adopted in our simulations. 

339 Figure 2. Relative importance of exospheric and surface sources of lunar pickup ions, (a) 

340 Flux range of ions from the exosphere computed at a dayside altitude of 100 km over two 

341 selenocentric locations, the subsolar point (blue), and a location of maximum flux (red); 

342 (b) flux of ions from the surface due to solar wind impact (fi-om Elphic et al, 1991). 

343 Fluxes were computed for a Parker spiral IMF configuration of 8 nT, the solar wind 

344 speed of 400 km/s, and the density of 10 protons cm*^. 

345 Figure 3. Distribution of He"^ and Na"^ of exospheric origin on the lunar dayside (alt=100 

346 km): (left) omni-directional fltix; (right) flux that maps into a nadir pointing 2 ti detector, 
34'’ a proxy for the Kaguya IMA. Both surface and exospheric ions map onto a spacecraft 

348 north of the white dotted lines, but surface ions are excluded from this calculation. 

349 Figure 4. Exospheric Na'^/He^ ratio that would be locally measured at a dayside altitude 

350 of 100 km by (a) a 4jt detector; (b) a nadir pointing 2n detector. The velocity distribution 

35 1 function (VDF) of ions arriving at the detector, in normalized units, is shown as a 

352 function of the electric field direction over two locations at this altitude, (c) a location 

353 near equatorial dawn, and (d)a location near the subsolar point. 'The portion of velocity 

354 space not seen by the nadir pointing 2ji detector is shaded. Because the fraction of flux 

355 measured varies with species and location, (e) exo-ion spectra integrated along any 

356 meridian would lead to underestimation by Kaguya IMA (red line) of Na"^ relative to He"^ 

357 for this IMF orientation. 
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